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Motivated by recent work involving the graviton-graviton tree scattering amplitude, 
and its twin descriptions as the square of the Bel-Robinson tensor, Bnuaj3i and as the 
"current-current interaction" square of gravitational energy pseudo-tensors t a p, we find 
an exact tensor-square root equality B^ va p = d 2 t a p, for a combination of Einstein 
and Landau-Lifschitz t a g, in Riemann normal coordinates. In the process, we relate, 
on-shell, the usual superpotential basis for classifying pseudo-tensors with one spanned 
by polynomials in the curvature. 

1. Introduction 

This paper revisits and relates the three ancient subjects of our title in a novel way, suggested 
to us in a very different context, the construction of invariants of D=ll supergravity Q. The 
graviton-graviton tree level amplitude M, particularly in the present exclusively D=4 context, is 
remarkably simple, namely the square of the famous Bel-Robinson tensor [|[ , 

L ~ B %aP • O) 

To be sure, there is some work involved here: First, since M is generated by intermediate graviton 
exchange, it is nonlocal, but when multiplied by the Mandelstam variables stu, it reduces to a 
local invariant L. The latter must, by power counting and (abelian) gauge invariance to this lowest 
(ft^/^) 4 = {g^u — V^u) 4 order, be a scalar proportional to R^ ua p- [The external gravitons are all 
on-shell {R^ =0), so R^uap really means the Weyl tensor throughout.] Here the power of using 
a suitable basis to classify scalar or tensorial powers of curvature || first emerges. For example, 
there are seven algebraically independent R 4 scalar monomials in generic dimension, but these can 
be shown to reduce to just two in D=4, owing to identities such as R lia p 1 R l/a ^' y = ^^R^p^s valid 
only there. This basis is spanned by (E 2 ,P 2 ), where (E,P) are the D=4 Euler and Pontryagin 
invariants. Since graviton scattering must be maximally helicity conserving (as follows directly 
from the supersymmetrizability of Einstein theory), this singles out L ~ (E — P)(E + P); but 
another deep D=4 fact is that this is in turn just B 2 . 

*deser@brandeis . edu, jf rankli@aprilia. Stanford. edu, domenico . seminar a@lpt . ens . f r 

f Present address: Stanford University, PO Box 13236, Stanford, CA 94309, USA 



1 



The final objects in our title, the gravitational energy pseudo-tensors, enter in a very physical 
way as the currents that generate the scattering through the cubic vertex ~ kJ d A xhdhdh in the 
expansion of the Einstein action, which can be expressed as the coupling of the field h a p to a 
current t a/3 ~ K(dhdh) a @ . The amplitude can then be expressed in the usual current-current 
form M ~ J t Xa D\ ua pt a P , where D is the intermediate graviton's propagator. Although gauge 
invariance further requires inclusion of the quartic contact term k 2 J d 4 xhhdhdh, it is plausible 
that this characterization of M yields another version of L, 

L ~ {dl v t aP f (lb) 

upon keeping track of momentum dimensions. Hence the temptation to equate the "tensor square 
roots" of (la) and (lb), 

B^vasp = t a p . (2) 

Let us next recall some properties of t a p and earlier attempts to relate B and t. First, we 
note that gauge systems with spin >1 do not even possess gauge invariant stress tensors, but only 
integrated Poincare generators Q. While this local defect is irrelevant in the absence of gravity, 
it becomes critical when the stress tensor is to be its source, not least for (self-coupled) general 
relativity itself, as witness the history of its countless energy pseudotensors (see e.g., ||). The 
problem is of course the clash between the second derivative order required to build the (linearized) 
in- or (generically) co- variant curvatures and the first derivative order building blocks of the t a g. 
Indeed, the Bel-Robinson tensor arose as an attempt at constructing a "more covariant" tensorial 
quantity modelled on the Maxwell stress tensor. However, the price is high: B has too many 
derivatives and indices, it is only covariantly conserved and it is not a physical current ||. Clearly 
any putative relation like (2), even at linearized level, must be a non-invariant one. Nevertheless, 
there exists a textbook exercise |7j proposing the equality B^ va p = d^ v t a p + Sp Va p in Riemann 
normal coordinates (RNC) for a particular (Einstein) t a p, but with a rather mysterious remainder. 
In achieving (2), we will have to consider all possible pseudo-tensors; as we have mentioned, there 
is an infinity of them and being identically ordinarily (rather than covariantly) conserved on shell, 
they can all be expressed as superpotentials there. This is also the case for B: on shell, it reduces 
to the identically (covariantly) conserved Sachs tensor [^]. In establishing (2), in the particular 
RNC frame, we will in fact translate the superpotential "basis" for the t a p and for B into one for 
four-index curvature quadratics. This will enable us not only to verify that it can be done there, 
but also how uniquely. 

Despite our arguments in favor of the existence of an exact connection (2) for some t a g, its 
validity is far from obvious; B is totally symmetric and traceless, while the d^ u t a p do not even 
display (a/3) symmetry in general, let alone the other invariances of B. It is therefore amusing that 
this correspondence, even if highly gauge variant, can be established at all! 

2. Ingredients 

A. Bel-Robinson. As mentioned, B first appeared in the endless search for a covariant version of 
gravitational energy density; the analogy with the Maxwell stress tensor T^ v = F m F v a +*F^ a *F U a 
(*F is the dual field strength) is striking in either of the two equivalent expressions, 

B^ii/aP a Rpuaf3 Bf„ a Rp Ua p , (3) 

Bfj,i/a(3 — B-^fj, aBpuaf) ~\~ R p pRpvaa ~29iivRa Rftpor ■ (4) 
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Here the dual curvature is *W ,V \ a = |e M " a@R Act- The interplay between these two expressions 
underlies the various special properties that B enjoys in four dimensions, due to on-shell identities 
that (like the dual itself) are only valid there. These include, besides the (fi <-» v) and (a <-> 0) 
symmetries, two further ones: (a) (fJ,v) (a/3) symmetry. This is transparent in (^); requiring it 
in (||) demands that 

9iivR(£ aT Rppor = 9apRff aT R-vpar 1 (5) 
i.e., that R a ^ T R^ par be a pure trace: 

R-a" 77 RppcJT — \gaf3R fipCrT Rfipar = . (6) 

Thus, this well-known D=4 identity is encoded in B. (b) (// *-* u) symmetry: "dualizing" B with 
e "/rpu j n ^ obviously annihilates it by the "** = —1" property. In terms of @), this then implies 
the more general 4-index relation 

^RfJ-a ^ Ryj3pcj R p fiRpucra R^v pRpaafl gidpvdaft SpfSdiya^RxpcrrR ^ — (7) 

of which (||) is the (a/3) trace; indeed these identities "explain" why B is necessarily 4-index, rather 
than 2-index like t a p. [For D>4 the identities (6), (7) cease to hold and a three-parameter family of 
conserved U-like generalization (each with partial properties) can be constructed; for more about 
B's and the curvature quartics BB in D>4, see pL] 

B. Riemann Normal Coordinates (RNC). We recall that at any one point in a Riemann space, 
coordinate invariance can be exploited in order to simultaneously (i) rotate the metric to Minkowski 
form (local flatness), (ii) annihilate all affinities (free fall) and (iii) remove 80 of the 100 d^ u g a p 
by use of the 80 independent <9^ 7 £<5 components of the gauge functions £5 in the curvature, leaving 
only 10 d 2 g combinations to represent the 10 Weyl curvatures R^vafi of interest on shell. In addition 
to (i) and (ii), 

9pu = 7]pu, g» u = rf u , g„u,a = 0, (8) 
the choices (iii) defining RNC are summarized as follows: 

— 3 gpv,a/3 = Rpavfi + Rpfiva , ~ 3 ^pu,/3 = R° fiw/3 + R" vp(3 ■ (9) 

Note that raising and lowering indices "passes through derivatives", here denoted by commas. 
On-shell, this means in addition, that all traces of (9) vanish: 

sTg^afi = o , = o , = o = r^ jQ . (10) 

C. Pseudo-tensors. As we have stated, these can be parametrized either by all possible independent 
superpotentials (of second derivative order), or more usefully in RNC by the appropriate curvature 
basis. Before proceeding to the general case, however, we write directly the two specific t a p that 
will enter in (2), namely the non-symmetric Einstein E% and symmetric Landau-Lifschitz L al3 . 
This will also give an idea of how RNC simplifies matters. We only keep those terms in each that 
will not vanish on-shell because of (10) after taking the two further, d^ u , derivatives: 

^ = -2r^r^ + ^ T r; A (n a ) 
v# = -r^ A + r^rt - (r^rj A + r^ A ) + g ^T x aT n T . (lib) 
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In writing these expressions, we have used the economy of notation permitted by RNC: first, since 
both E and L are special in being bilinear in T's (i.e., not involving terms like gd 2 g), then both of 
the d 2 u must act on these T's and not on any undifferentiated metrics; we have therefore set all the 
latter to their Minkowski values, so that summation and moving indices in (11) is to be understood 
in light of (8). Note also that, unlike L a @, Eg (here expressible as E a @) is not yet (a/3)-symmetric; 
this will, however, come about after differentiation. 

D. Quadratic Curvature Basis. All 4-index curvature bilinears are algebraically equivalent to the 
double contraction 

where we have indicated the symmetries of Q; there are no further ones, either within a pair or 
under index exchange between pairs. The three basis members we need are representable by 

— Qafj,uf3 ^jnuafi — Qaf3fiv Z[iva(3 = Qa/if3u (+ l/ / i ) > (13a) 

where (+z^) means that each quantity is to include the (^/^-interchanged form in its defini- 
tion. These three linearly independent quantities are a subset of the complete (6-member) non- 
symmetrized basis and agree with the enumerations in Q. It is convenient, in addition, to define 
a single trace object, which we take to be 

T^i/af) = — \ 9/iu Qaaaf3 = dfivSafl R-XajS ■ (13b) 

All the (X, Y,Z,T) are uniformly labelled by the indices (/j,ua(3) in that order; we will avoid index 
proliferation below by (usually) leaving them off altogether. In terms of this basis, B^^p is the 
combination 

B = Z + 3T . (14) 

It is also convenient to define the tensor S^ va p = R^P a R„/3 pcr + R^Rpppa + \9iiu9afiR\bcd which 
is expressible as 

S = 2B-2X -12T . (15) 

It also follows by (7) that 

S = 2Y -2B-12T + 6(T^ + ) ■ (16) 



3. The Relation 

In RNC, any d 2 u t a p will be a sum of products ~ (d 2 g d 2 g) a p pu , with manifest {(iv), but 
not necessarily (a/3), symmetry. There is thus a small class of possible terms; all are expressible 
in terms of the basis (X, Y, Z, T) of (13). We start with the Einstein contribution (11a); using the 
prescription (9), we obtain 

\dl u E a p = \X + lZ + T . (17) 
In terms of B and S this can be rewritten as 

\df Ll/ E al3 = ^(B -\S), (18) 

the original relation proposed in Q. The d 2 w L a g of (lib) contains more terms but is equally 
straightforward. We find for it 

dl v L aP = -\X + lZ + 2T , (19) 
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or equivalent ly 

0jUa/J = +£(7B + ±S). (20) 
Adding (17) and (19) yields (14); hence the promised formula, 

Ba(3^u = <^L (L a p + | E a f^j . (21) 

In contrast to (18), there is no remainder here. Incidentally, although (21) has been derived for 
D=4, where there is a unique B, it can be extended to any D using the B of (4). 

There are two immediate questions about (21): How unique is our result, and if it is, why 
just this combination of t a g? We have no wisdom as to the latter: the invariant amplitude provides 
no hints about why precisely (L a p + \ E a p) requires no remainder in this frame. To quantify the 
uniqueness aspect, we will write down all superpotentials in RNC, since these cover all t a p. There 
are both symmetric and non-symmetric A Q/ g (relevant because we used E a p ^ Ep a ). Symmetric 
ones are of the form S a @ = d\ a H aXl3a , where H has the algebraic symmetries of the Riemann 
tensor; this ensures both identical conservation and symmetry of S a ^ . Nonsymmetric ones will be 
conserved on only one index, A a ^ = d\H aX @ where H aX @ is only antisymmetric in (aX). Since 
we are concerned in RNC with terms of the form d 2 g d 2 g and in particular with d 2 „S al3 , H aX ^ a 
depends only on the metric, while H aX @ will be ~ gdg. As shown in the Appendix, there just are 
enough independent superpotentials to span the complete (X, Y, Z, T) basis and thus to express B 
uniquely. Hence, there is just one effective t a p - in RNC - that fulfills our relation, though t a( g's 
that are intrinsically different in an arbitrary frame may degenerate to a single one in RNC. 

4. Summary 

Our modest result is that on shell, at the origin of the RNC frame, there is an exact local 
equality between the Bel-Robinson tensor and the double gradient of a particular energy pseu- 
dotensor. We have conjectured this relation to be the "tensor square root" of a more physical one 
in tree level graviton-graviton scattering, whose amplitude is simultaneously proportional to the 
square of B^ ua p and (essentially, if not quite gauge invariantly) to that of a pseudotensor. Despite 
all this fine print, one cannot help but wonder if there is more to be learned from (21). 
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Appendix: Superpotentials 

We enumerate here the independent superpotentials in RNC. Our purpose is to show that 
they in fact constitute a basis equivalent to the (X, Y, Z, T) , in terms of which we know that our 
class of tensors is expandable. Consider first the symmetric ones; we are actually concerned with 

AJU* = 9lu S a p = d% Xa ff**>{g) . (A.l) 
As explained in text, H contains no derivatives. It is then of the form 

H aX(Sa = {g a pg Xa - g a *gf]x)K(g) (A.2) 

where K is a "scalar" like 1 or sf—g- The four derivatives will then all fall on the (gg — gg) term, or 
all on K(g) to form the surviving d 2 d 2 g combinations, but not two on each: a d 2 K(g) is excluded 
because it would necessarily be proportional to g Xr7 d 2 g\ a which vanishes on shell, irrespective of 
the d 2 indices. When all four derivatives fall on the (gg — gg), the undifferentiated K is a constant, 
K(g = rj^), and it is easy to see that this gives rise to the contribution 

9A^ = (-5X + 2Y + 4Z)K(ri) . (A.3) 

When all derivatives act on K alone, we have the usual transverse projector in (a/3), 

A^ = (r^D - dl p ) d% (g x °g Xa ) (A.4) 

where we have kept the relevant, quadratic in metric, part of K; two derivatives are to be distributed 
on each metric, neglecting terms such as Og Xo d 2 u g\ a that vanish on shell. The significant point for 
us is the presence of a (unique) trace term here, i.e., that A| includes a part 

VafSg^garM ~ ? , (A-5) 

using (9), the identity (6) and the definition (13b). As we shall next see, there are enough inde- 
pendent antisymmetric contributions to span the remaining (non-trace) basis members (X, Y, Z) 
between them. The three possible forms are given by 

A a\f3 = g\agal3,a ~ g a ag\l3,a ^-aXfS = (gXa,a ~ gaa,x)g/3a ^ a X/3 = gaag/3cr,X ~ gXag/3a,a ■ (A.6) 

In principle there can again be coefficients K(g) here, but in fact they will not contribute: We are 
interested in expressions ~ d^ va (gdg K). Only the d 3 (gdg) K(rj) part fails to vanish: But as in the 
symmetric case, d 2 K(g) ~ d 2 (g Xcr g\ a ) always vanishes on shell. Hence we have the three possible 
terms A\ = d^A^ with the A { of (A.6). The result is as follows 

9A\ = -5X + 2Y + 4Z = 9A{ , 9A 2 A = 3X - 6Z , 9A 3 A = -4X + Y + 5Z . (A.7) 

Being linearly independent, these A* are equivalent to (X,Y,Z); they contain no trace, which is 
instead carried by A|, as explained above. Thus, in RNC the set of four possible superpotentials, 
like the (XYZT), span all (double gradients of) the pseudotensors t a p as well as B. 
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